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ABSTRACT 

Binary populations in young star clusters show multiplicity fractions both lower and up to twice as high as those observed in the 
Galactic field. We follow the evolution of a population of binary stars in dense and loose star clusters starting with an invariant 
initial binary population and a formal multiplicity fraction of unity, and demonstrate that these models can explain the observed 
binary properties in Taurus, p Ophiuchus, Chamaeleon, Orion, IC 348, Upper Scorpius A, Praesepe, and the Pleiades. The model 
needs to consider solely different birth densities for these regions. The evolved theoretical orbital-parameter distributions are highly 
probable parent distributions for the observed ones. We constrain the birth conditions (stellar mass, M cc i, and half-mass radius, rj,) 
for the derived progenitors of the star clusters and the overall present-day binary fractions allowed by the present model. The results 
compare very well with properties of molecular cloud clumps on the verge of star formation. Combining these with previously and 
independently obtained constraints on the birth densities of globular clusters, we identify a weak stellar mass - half-mass radius 
correlation for cluster-forming cloud clumps, n,/pc oc (Me C i/Mo) ' l3±004 . The ability of the model to reproduce the binary properties 
in all the investigated young objects, covering present-day densities from 1-10 stars pc~ 3 (Taurus) to 2 x 10 4 stars pc~ 3 (Orion), 
suggests that environment-dependent dynamical evolution plays an important role in shaping the present-day properties of binary 
populations in star clusters, and that the initial binary properties may not vary dramatically between different environments. 

Key words, binaries: general - star clusters: general - globular clusters: general - open clusters and associations: general - meth- 
ods:numerical 



1. Introduction 

It has been observationally established that the properties of bi- 
nary populations in star clusters depend on environment: the 
frequency of binaries is lower in denser and older populations. 
For instance, the Taurus star-formation region, which is about 
1 Myr old, has a very low density (1-10 sta rs pc~ 3 ) and a bi- 
nary population th a t is about twice the size dGhez et al.l 1993; 
i Leinert et ail 1 19931; Isimon et al.|[l99l iKohler & Leinertlll998l 
■ iKraus & Hillenbrand! l2009h of those for solar-type main se- 
quence stars in t he Galactic field (the canoni cal sample used here 
for comparison, Duquennov & M avorll99li) . On the other hand, 
the similarly aged Orion Nebula Cluster has a much higher den- 
sity (2 x 10 4 star s pc~ 3 ) and its binary population is smaller than 
that in the field dReipurth et al.ll2007l) . The density of IC 348 is 
comparable to that of NGC 2024, but the former is 3 - 5 times 
older. IC 348 has a smaller bina ry fraction than NGC 2024 
dDuchendl 19991: iLevine et"ai]|2000|) . 

Similar behaviour is rep orted for older open (Sollima et al.l 

120101) and globular clusters (ISollima et al.ll2007l) of the Galaxy, 
inferred from anti-correlations between the binary proportion 
and cluster age, as well as between binary proportion and cluster 
luminosity (i.e. mass). 
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Theoretical approaches to explaining these differences can 
be categorized into two major groups. One proposes that the 
formation of binaries dependends on the environment. In this 
case the prevalent conditions in the star-for ming cloud can de- 
termine the outcome of bi nary star formation (iDurisen & Sterzikl 
119941: ISterzik et al.ll2003l) . The second group of scenarios pro- 
pose environment-dependent dynamical evolution, in which the 
properties of binary populations seen in star clusters cannot be 
assumed to be primordial, i.e. they are not directly related to the 
birth properties of binaries. In these scenarios, the characteristics 
are modified by interactions between the stars and binaries. Even 
if the clusters are still very young (< 1 Myr), the binary prop- 
erties can strongly deviate from the birth conditions produced 
in the star-formation process since the binary proportion is re- 
duced rapidly and orbital paramet ers are altered efficiently over 
only a few initial crossing-times (|Kroupall995al : iDuchenef 1999; 
iFregeau et al.ll2009HParker et al.ll2009l) . 

In this paper, we investigate whether dynamical evolution is 
able to explain observations of the multiplicity properties in very 
different environments, starting wi th an assumed u niversal initial 
binary population in star clusters (lKroupalll995bl) . In particular, 
we apply the concept of inverse dynamical population synthesis 
to constrain the properties of star clusters at their birth. 

In Sect. |21 we present the model and the procedure to analyse 
the eight investigated regions we describe in Sect. [3] In Sect. [4] 
we present the outcome of our analysis and Sect. [5] summarizes 
and discusses our results. 
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2. Binaries in star clusters 

2.1. Stimulated evolution, dynamical equivalence, and the 
age-density degeneracy 

A stellar binary population in a star cluster evolves as binaries in- 
teract gravitationally with other stars or binaries. As energy and 
momentum between these centre-of-mass systems is exchanged, 
binary orbital parameters are altered or, if the transferre d amount 
of energy is sufficient, a binary dissolves. Generally during en- 
counters, a soft binary, which has an internal binding energy 
that is lower than the kinetic energy of a typical centre-of-mass 
system in the cluster, becomes on average less bound (softens), 
while a hard binary with a higher binding energy be comes more 
strongly bound (hardens, lHeggielfl975l: lHillslll975l) . This stim- 
ulated evolution of binaries depletes the binary population. The 
process is most efficient in wide low-energy systems, which are 
those most likely to be disrupted. 

The resultant change in the orbital parameters for popula- 
tions of binaries evolving in star cluste r s drive n by stimulated 
evolution is quantified by iMarks et all (1201 ll) by means of a 
comprehensive set of A^-body computations of binary-rich star 
clusters. In these models, all stars are initially members of a 
binary. The component masses of late-type stars are selected 
randomly from the canonical stellar initial mass-function (IMF, 
iKroupaO OOl). The IMF is conveniently described as a two-part 
power-law in the stellar regime, 



P fid =10 3 °M 



£(m) oc m a ' where 



ai = 1.3 0.08 < m/M G < 0.5, 
a 2 = 2.3 0.5 < m/M Q < 150, 



(1) 



with an average stellar mass m » 0.4 where dN = %(m)dm is 
the number of stars in the mass interval [m, m + dm). Random 
pairing results in late-type binaries having a flat initial mass- 
ratio distribution which is affected only in a minor way thr ough 
pre-mai n sequence e i genev olution (fig. 17 in lKroupal d2008l) and 
fig. 2 in lMarks et alj (1201 ll) ). The mass ratio distribution of late- 
type binaries remains approximately flat after dynamical pro- 
cessing. The orbital periods of the binaries are selected from an 
initial period distribution function that increas es with increas- 
ing period and flattens at the longest periods (lKroupal [T995b. 
Fig. Q]). This distribution agree fl with pre-main sequence and 
pro to- stellar multiplicity data (jkroupa & Petr-Gotzensl 120111: 
Marks et al. 201 1). The period distribut ion is assumed to be the 
universal outcome of binary formation (lKroupall201 ll) . 

Allowing this invariant initial binary population to evolve 
for 5 Myr i n clu sters of different initial masses and sizes, 
Mark s et ail (1201 ll) find that for their set of models (of stellar 
mass in the computed clusters M ec \ < 1O 35 M and half-mass 
radii r n < 0.8 pc), the evolution is driven by the initial cluster 
density, p ec i = 3M ec i/87rr 3 , i.e. models of the same p ec i evolve 
their population in the same way. They are dynamically equiva- 
len t (over the 5 Myr int egration time). 

IMarks et"ail (1201 ll) quantify a stellar dynamical operator, 
Q dy ™'' rh (f), which transforms an initial orbital-parameter distri- 
bution, (e.g., x =semi-major axis or x 

F.M,. c |,r h 



evolved distribution, <D r 
evolution 



mass-ratio), into an 
(f), after some time, t, of stimulated 

(2) 



1 Note that this initial distribution is not a fit to any observed data, but 
was identified by inverse dynamical population synthesis (Sect. l2~2l . It 
is however based on an initial guess of a flat distribution that matches 
the pre-main sequence constraints. 
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Fig. 1. The Q-operator in action (eq. [2] thick solid S-shaped 
curve, schematic, right ordinate). Its exact shape depends on 
the cluster volume-density within the half-mass radius, p ec i, and 
time, t, for stimulated evolution. Multip lying the resp ective oper- 
ator by the initial binary population (Krourja 1995b, thin dotted 
line) for a cluster with p ec i = 1O 3 M pc -3 , and age, t = 1,3, 
and 5 Myr, respectively, results in the three semi-major axis dis- 
tributions (solid histograms from top to bottom). Most of the 
evolution has already taken place in the first Myr, where the 
binary fraction decreases from 100 to 72 percent in this par- 
ticular cluster. The grey thick dashed line on top of the lowest 
black solid line is the resulting semi-major axis distribution for 
a Peci = 1O 3,5 M pc~ 3 cluster after 1 Myr of stimulated evolu- 
tion. This shows that two clusters with densities p ec i(l) < Peci(2) 
can evolve into the same resulting distribution when we assume 
f(2) < f(l) is appropriately chosen (owing to the age-density de- 
generacy). That is, the p ec i(l) cluster at age f(l) is dynamically 
equivalent to the p ec i(2) cluster at age f(2). 



Eq. (f2]i allows the efficient calculation of the evolution of a bi- 
nary population in a star cluster of a given initial density without 
the need for expensive individual A^-body computations. Fig. [TJ 
depicts how the operator works. The evolution of the binary 
populatio n has largely ended by 5 Myr, the time for which the 
models in Marks et al] d201 ll) are computed. We note that there 
is an age-density degeneracy, i.e. a high density cluster of age 
f(l) may have the same binary population as a somewhat lower- 
density cluster of age f(2) > f(l) since stimulated evolution in 
the latter cluster is less efficient. 

In this paper, the considered orbital parameter distribution is 
always the distribution of semi-major axes, log 10 a, of the con- 
sidered region. For a range of primary masses around m\, it is 
defined as the distribution of binary-fractions /bin = Nb/N cms as 
a function of log 10 a 



^M ec] ,r h , . t//bin(log 10 fl) 1 dN b (\og m a) 

O. (mi) = = 

l °Sio" d\og m a N cms c/log, n a 



(3) 



'10 



where iVj and Af cms are the number of binaries and centre-of- 
mass systems (singles+binaries), respectively. The total binary 
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fraction equals the area below the distribution, i.e. the integral 
over the semi-major axis distribution is 



/bin 0*1 ) 



J lo gio " 



(mi) d\og m a 



(4) 



For a description on the method and extraction of distributio n 
functions, the reader is referred to Marks, Kroupa, & Oh (2 01 lh . 



2.2. Inverse dynamical population synthesis 

In this paper, the inverse operator, (Q^ 1 ' a )~ 1 ( f X is applied to 
the observed binary distributions in clusters of age t in order to 
see whether stimulated evolution in a cluster of initial density 
Peci can turn the universal initial distribution into a shape which 
serves as a parent function of the observational data. 

This inverse dynamical population synthesis (IDPS) relies 
on assuming an invariant birth binary star population and match- 
ing an observed binary population with a dynamical model from 
which the initial conditions of the population in terms of its den- 
sity are derived. Kroupa (1995a) introduced this concept and 
demonstrated that physically relevant solutions are obtainable. 
He postulated that two stellar dynamical operators, p £2^ l '' h and 

q Q!^' Vh , exist. The first one ( p Q.^' rh ) evolves the period distri- 
bution function of a binary population through stellar-dynamical 
encounters, while the other one CfiJ' 11 ) evolves the mass- 
ratio distribution. As both distribution functions are not corre- 
lated for periods > 10 3 days, solutions for p nf od,n " and «Qf 

r J dyn dyn 

can only b e obtain e d if the se operators correspond to physical 
processes. iKroupal d!995al) not only discovered the existence 
of these operators, but also demonstrated that they are equal, 
i.e. the same process (stimulated evolution) 



dyn 



= m 



M cc ,,r h 
dyn 



'dyn 



ex- 



shapes these distributions. Thus, a single operator, Q 
ists to turn both the initial period and the initial mass-ratio dis- 
tribution into the observed shape simultaneously. It is given by 
a star cluster with M ec i = 128M and r n = 0.8 pc, which cor- 
responds to the typical embedded cluster from which the major- 
ity of Galactic fi eld binaries originate (cf . also Marks & Kroupa 
120111) . [Kroupa (200l and iMarks et all (|2011|) formalized the 



mathematical concept of 0^' rh , which is here used for the first 
time to derive the initial conditions of observed star clusters. 



2.3. Modelling observed binary populations 

We construct the observationally inferred semi-major axis distri- 
bution, (D° bs , from N orbits distributed into «bin bins, and the 

log 10 a 

evolved theoretical distribution, <&, ed, ' i (?), where t is chosen to 
correspond to the age of the investigated region for which the 
p ec i is known (eq. [2j. The best-fit solution for p ec i is then found 
by performing a^- 2 -test with «bi n ~ 1 degrees of freedom, i.e. the 



sum of the squared deviations between <S? bs and 4>, ed ' rb (f) is 

1 log 10 a log 10 a w 

minimized for p ec i 



"bin 



b,„ (of ed ' rh (f)-0° bs ) 

2/ \ V 1 V log 10 fl v ' log,o"7; 
X (Peel) = 

i=\ 



V logio a h 



mm 



(5) 



where AO? bs is the error in the observational data in bin ;, 

which is assumed to be V^V errors for all data sets. 

To estimate the likelohood that the resulting Of^ 'V?) is in- 

& logio a ' 

deed a parent function of O° b o s a , we calculate the probability, 



P(> xtm)' tnat a P oorer At occurs dPress et al.lll994l) . For in- 
stance, at a given v 2 . , P(> y 2 ■ ) - 0.3 means that if one re- 

' fc ^min' v ~ A -min / 

jects the null hypothesis (that "G^Jf '^(f) is a parent function for 



0° bs "), then there is a 30 percent chance that <bf< d ' rh (t) is a 

logio" r logio" 

parent function, even though it has been withdrawn. Confidence 
intervals around the best p ec i are defined by the condition P(> 
xtyj„) =10 percent, outside which the null hypothesis is falsely 
withdrawn with < 10 percent chance only. 

The,y 2 -test allows the identification of a typical birth density 
Peci for the considered region and via eq. (|4} the overall present- 
day binary fraction can be assessed using the best-fit theoretical 
distribution. 

While the dynamical operator acts on a population of bi- 
naries of all primary masses, the observations mostly probe 
masses of a restricted primary-mass range. As both the bi- 
nary properties and their distributions are mas s dependent 
dKouwenhoven et alj2009HMarks & Kroupal201 lh . this effect is 
important to consider. The analytical model allows us to extract 
only those binaries in the theoretical population that cover the 
sa me range of stellar ma sses as the observations (see sect. 2.4 
in lMarks & Kroupal201 lh . However the effect on the results is 
negligible as the observational uncertainties are much larger than 
the mass-dependent effects discussed here. 



3. The young cluster sample 

We now introduce the young cluster sample that we analyse with 
the model described in Sect. 12.31 For each cluster, we compile 
information about their present-day masses and sizes and their 
binary population. We emphasize that when discussing the dy- 
namical processing it is indispensible to distinguish between the 
dynamically active sub-clumps or individual clusters that always 
have radii of about 0.5 pc and the whole coeval stellar popula- 
tion that may extend over a few pc (see also the discussion in 
Sect.©. 

We ascertain from the literature whether the initial stellar 

mass distribution is compatible with the canonical IMF, which 

i 1 j 1 

is used as input into the A^-body models of Marks et al. (201 1). 
The theoretical distributions to compare with the observational 
data are evolved for the age of the respective region. Since these 
ages are often uncertain or an age spread is present, we use 
mostly upper age-limits, as indicated in Figs. [2]to|4] to construct 
the theoretical distribution. Owing to the age-density degeneracy 
(Fig. [TJ, this will result in lower limits to the inferred initial den- 
sities. In addition, when the ages are chosen to differ slighty the 
resulting change in the densities is smaller than the given con- 
fidence intervals, which are rather large given the observational 
uncertainties. Therefore, the results are considered reasonable 
estimates of the initial densities. 



3.1. Taurus 

The Taurus star-formation region has a d istance of 140 pc 
dKenvon et alJ 1 1994b IWichmann et alJ Il998h . which is nearby, 
and extends over a region of (100°) 2 on the sky. The age of its 
population is « 1 Myr. Six dis t inct cl usters of stars in Taurus 
were detected by Gom ez et al.l (Q993) to have projected radii 
of 0.5 - 1 pc with ss 15 stars or binaries (with a mass of 
about 6M for m » 0.4) each and a median separation be- 
tween the stellar centre-of-mass systems within each cluster of 
0.3 pc. The individual clusters are separated by a few pc, thus do 
not interact dynamically. Taurus contains more than 300 known 
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members and the stellar de nsity in each of its clusters is at 
1-10 stars pc~ 3 very low (lLuhmanetal.ll2009t) . The N-body 
models of iKroupa & Bouvien d2003f) show that Taurus Auriga 
like clusters are largely unevolved and in particular that the bi- 
nary properties should be similar to those at birth. 

Taurus was originally prop osed to be deficient in low- 
mass stars (iBriceno et al.ll2002l). but many brown dwarfs hav e 
since been discove red (ILuhmanl l2004b|: iGuieu et aD |2006). 
IKroupa et all (120031) and lThies & Kroupal 00081) show that there 
is no need to invoke a non-canonical IMF in Taurus. 

Am ong the 43 observed primary stars in Taurus, iGhez et all 
(1993) find 16 binaries with angular separations between 16 - 
252 AU. The binar y fraction i s thus 37 ± 9 percent over this 
separation range. Simon et al. (1995) find 22 binaries and 4 
triples among 47 systems. The binary frequency in the range 
3 - 1400 AU is at leas t 1.6 + 0.3 times the value of the canonical 
Galactic field sample. iRraus & Hillenbrand! d2009l) find 27 wide 
binaries in Taurus in the separation range 500 - 5000 AU, which 
is consistent with a log-flat distribution but inconsistent with the 
Galactic field log-no rma l distribution. Here the co mbined data of 
iLeinert et al.1 (H993) and Koh ler & Leinerj (1 19981) are used, who 
identified in total 74 binaries or multiples among 174 systems. 
The binary fraction of 42.5 + 4.9 percent is larger by a factor of 
1 .93 + 0.26 than that of solar-type main sequence stars over the 
accessible separation range. 

3.2. p Ophiuchus 

The p Ophiuchi cloud is a complex of dark nebulae and molec- 
ular clouds with ongoing star forma tion located at about » 
120 - 160 pc dista nce from the Sun dChinilll98lHKnude & Hogl 
1998t lMotteeTaIi ri998. 140 pc is adopted). Its densest part is 
the L1688 dark cloud, which is itself structured into the six 
major clumps Oph-A to Oph-F, each containing a few tens of 
solar mass in young stel lar objects (YSOs) and each having 
a diam eter of ss 0.3 pc ( Stamate llos et alJ 120071) . lLoren et al] 
(1 19901) use a total mass (YSOs + gas) for each clump of 
=s 2OOM . The a ges of the stars are inferred to he in the 
range 0.3 - 1 Mvr dWilking et al.lfl989l: iGreene & Meverll 19951 
iLuhman & Riekelll999l) . iBontemps et al] d200ll) find, after cor- 
recting for an assumed binary fraction of 75 percent, a stellar 
mass-function for young stellar objects (YSOs) that is consis- 
tent with the slope of the canonical stella r IMF for low masses 
(0.055 - 0.55M o ) and agrees with earlier dComeron et al.lll993l 
IStrom et al.ll 19951 iMotte et alj|l998t iJohnstone et al.ll2000l) and 



later dTachihara et al.ll2002HStanke et al.ll2006l) studies. A some- 
what steeper mass-function than the canonical I MF is observed 
for h igher stellar masses (012 - 1 .7 in eq. Q] Bontemns et al] 
120011) . The same authors estimate that there is 92M of mass 
in YSOs in L1688 and given an average radius of 0.4 pc for the 
whole region they find a stellar mass density of 34OM pc~ 3 . The 
extents of the clumps are about 0. 1 — 0.2 pc. The total gas mass 
of the re gion is between 550M^ dWilking & Ladal 1 19831) and 
1500M G dLiseau et all 1 19951) . iRatzka et al.1 d2005l) state a value 
of 600M o for the total mass within a region of 1 x 2 pc. 

The YSO content of L1688 has been intensively surveyed in 
the past two decades. IGhez et all d!993l) observed a binary fre- 
quency for T Tauri stars of 29 ± 11 percent (sample size of 21 
stars only) in the linear separation range 1 6 - 252 AU, which 
is significantly higher than that in the field. Simon et al. I (fl99l 
find 10 binaries, 2 triples, and 1 quadruple among 35 T Tauri 
primaries in p Oph in the separation range 3 - 1400 AU, i.e. a 
binary fraction of 37 ± 10 percent. Since they do not perform 
a completeness correction, this is at least 1 . 1 times larger than 



that in the field o ver the same range . iDuchenel d 19991) correct 
IGhez et all dl993l) 's and ISimon et al] dl995l) 's data for incom- 
pleteness and re-evaluat ed the companion-sta r fraction of p Oph 
to be 26 ± 5 percent. iBarsonv et all d2003l) find 24 + 11 per- 
cent withi n 14 - 154 AU amo ng 80 primaries in their restricted 
sample. Ibuchene et al] {2004) observe 41 YSOs in a deep near- 
infrared survey, finding 14 companions to 11 primaries, i.e. a 
companion star fraction of 34 + 7 percent. 

In this work, the data bv lRatzkaet al] d2005l) are used, hav- 
ing observed the to-date largest sample of YSOs in p Oph for 
their multiplicity properties. They find 112 single, 43 binaries, 
and 3 triples among their 158-object-sized sample, resulting in a 
binary fraction of 29. 1 ± 4.3 percent in the separation range from 
18 to 900 AU. This is slightly larger than the corresponding field 
value (23.5 ± 4.8 percent). Their semi-maj or axis dist r ibution 
is indistingu i shable fro m those published in | Ghez et al .1 d 19931) . 
ISimon et all d 19951) . and lBarsonv et al] d2003l see Fig.l2l lower 
left panel). 

3.3. Chamaeleon 

The star-formati on region Cha maeleo n is 160 - 170 p c away 
from the Sun dWhittet et al] 119971: IWichmann etafl 119981: 
iBertout et al] 119991) and has a median age of about 2 Myr 
dLuhmanl2004al) . Cha I and Cha II are the main cloud complexes 
in Chamaeleon where the majority of stars are being formed. The 
IMF is compatible wit h that in the field as well as IC 348 an d the 
Orion Nebula Cluster dComeron et al.l2000HLuhmanl2008l) . The 
total mass in the dark clouds of Chamaeleon is about 5000M o , 
Cha I containing 237 kn own members, and 50 members having 
been identified in Cha Il dLuhm an 2008). The total stellar mass in 
Cha I is about 100M H , most of it being contained within a radius 
of 1° (Luhman 2007), corresponding to 2.88 pc at the distance 
of Chamaeleon. 

Among 77 observed primaries, |Kohlerj(|2P01) find after cor- 
rection for background stars 10.3 binaries and 0.5 triples, respec- 
tively. The binary fraction in the separation range 10-480 AU is 
thus 14 percent, about nine percent smaller than the 23 percent 
determined for solar type stars in the Galactic field. Interestingly, 
this deficit occurs only in their two largest separation bins, while 
the smaller separation bins neatly agree with the observed field 
separation distribution. Thi s defici t of wide binaries is also re- 
ported by iBrandeker et al.1 (|2006). The most recent study of 
lLafreniere et al.l d2008l) . which finds 30 binaries and 7 triples 
among 126 primaries, does not however suggest that there is 
a strong depletion for the widest binaries. These authors pro- 
vide projected angular separations for their observed binaries 
and triples, but they do neither sho w nor discuss a separation 
distribution. Lafre niere et al.l d2008l) , s data is used here and the 
distribution agrees also at the largest separations with the field 
distribution. 



3.4. Orion Nebula Cluster 

The Orion Nebula Cluster (ONC) is a young cluster (< 1 Myr) 
that is pa rt of the Orion complex at a distance of about 414 - 
343 pc dMenten etal]l2007b iKraus et al] [20071) from the Sun. 
With its high central density, it is probably the denses t nearby 
star-formation region. Hillenbrand & Hartrriann] (11998b review 
earlier work on the ONC and determine its main properties. They 
recover 3500 centre-of-mass systems in total out to a radius of 
2.5 pc. They calculate that if the cluster were in virial equilib- 
rium a total mass of * 45OOM would be needed within 2 pc, 
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which differs from the observed stellar mass within this radius of 
1800M o . In the core of the ONC (size 0.16-0.21 pc), the stellar 
density is about 2 x 10 4 stars pc~ 3 . Stellar masses range from 0.1 
to 5OM with an average mass of O.8M with most massive stars 
being found in the centre. The core contains the Trapezium sys- 
tem consisting of ( at least) 6 stars with a total mass of 88.4M 
dHillenbrandl 1 1 9971) . Brown d warfs have since been detected 
in the centre of the cluster ([H illenbrand & Carpenter 2000; 
iLuhman et aill2000t Eucas et al.ll200U IMuench et al.ll2002l), s o 
that the mean mass may decreases to 0.5M^ dScallv et alJl2005h . 
in good agreement with the canonical IMF. 

Existing dynamical TY-body models suggest that the ONC 
must h ave been far m ore compact initially to appear in its present 
shape. [Kroupal d2000b computed models that start in virial equi- 
librium and are either dynamically hot (expanding) or cold (col- 
lapsing), with the same realistic initial binary population used 
here and expulsion of residual gas. He searched which initial 
conditions reproduce the observed velocity dispersion and ob- 
served number of stars in the core simultaneously. The mod- 
els initially in virial equilibrium were basically excluded, owing 
to the large number of stars needed to match the observational 
properties. Models starting with 3500 to 10000 stars initially 
(corresponding to 1400 to 4OOOM form = O.4M ) that collapse 
and expand allow for ONC progenitor s with half-mass ra dii of 
0.5 - 1.2 pc and < 0.5 pc, respectively. Scallv et al. (2005) com- 
pute very similar models to the iKroup a (2000) ones but with- 
out gas and binaries. Their computations agree in their main re- 
sults, although virial equilibrium models most closely reproduce 
the observed density profile. In this case, the initial density may 
have been one to two orders of mag nitude higher than that today. 
Pflamm- Altenburg & Kr oupa (2006) investigate the decay of the 
ONC core. Models with initially 40 massive stars (> 6OOM ), 
which is the expected number if stars in the ONC were selected 
from the canonical IMF at birth, in a region of 0.025 pc only end 
up with about the observed number of massive objects (10 stars 
above 5M ) and display a Trapezium structure. 

The binary population of the ONC has been found 
be about half as large as that in Taurus ( Prosser et al] 
1994t iPadgett et alJ Il997t iPetr et al.1 [19981: ISimon et al.lll99' 



to 



McCaughreanl 120011) . 
a higher binary free 
( Preibisch 



binary treq 
et al.lll999b [' 



The high-mass population exhibits 
uency than the l ow-mass primarie s 
tl et al.1 l200l IKohler etail |2006l) . 
For the ONC, we use the study of iReipurth et al.l (|2007), who 
find 75 binaries and 3 triples among 781 ONC members. In the 
limited range from 67.5 to 675 AU, this corresponds to a binary 
fraction of 8.8 + 1.1 percent after correcting for contamination. 
The field binary- fraction for solar- t ype stars in this range is about 
1 .5 times larger. iTobin et all (120091) add short period binaries to 
previous studies by using spectroscopic data. They infer a bi- 
nary fraction of 1 1 .5 percent in their sample up to a period of 
roughly 4000 days, as opposed to » 17.5 percent in the field. 
Unfortunately, however, they cannot assess their completeness 
or even report reliable orbital periods. 



3.5. IC 348 



IC 348 has a total stellar mass of 1 1OM dPreibisch et al.ll2003l) . 
a half-mass radius of 0.47 pc and an overall stellar density within 
this radius of 22OM pc 3 (lLada & Ladal fl99l . The IMF in 
IC348 is Salpeter-Uk e above O.7M and flattens below that value 
(IMuench et alJl2003[). which agrees well with the field-sta r (i.e. 
canonical) IMF dLuhman et alj|1998tlPreibisch et alj|2003l) . The 
cluster has a mean age of 2 Myr, but an age spread of * 3 Myr, 



and is located at a di stance of 320 pc from the Sun (lHerbigl l998: 

IMuench et al "1l2003l) 

iDuchene et al. (1999)'s multiplicity data is used in the 
present work. They investigated the multiplicity properties in a 
near-infrared survey of 66 low-mass members of IC 348, find- 
ing 12 binaries in the separation range « 40 to 3200 AU. The 
binary fraction of 19 + 5 percent agrees within the errors with 
23 ± 3 percent in the field over the same range of separations. 



3.6. USco-A 

Upper Scorpius (USco) is the youngest (« 5 Myr) region in the 
Scorpius-Centaurus OB association located at about 145 pc dis- 
tance, which contains USco-A. USco-A's stellar IMF between 
0.1 and 2OM is consistent with the field (i.e. canonical) stellar 
mass func tion and the total stel lar mass is estimated to be about 
2O6OM dPreibisch et al.i r2002). The same authors estimate the 
total initial mass (gas+stars) to be ss 80000M o and the initial 
diameter to be 15 pc given the p resent-day velocity dispersion. 

iBrandner & Koehleri fl998) carried out observations of 
68 TTauri stars in USco-A finding 21 binaries with angular sep- 
arations of between 15 and 4 35 AU, which is here used for the 
analysis. IKohler et alJ (|2000) probe a similar separation range 
(19 - 870 AU) but have a slightly smaller sample of 48 TTauri 
stars. A natural complement to this d ata is the wide binary sur- 
vey by iKraus & Hillenbrand! d2009h who identify 21 binaries 
betwe en 500 - 5000 AU. Here we use the data of IKraus et al.l 
(2008), who observe 82 systems finding 27 companion stars at 
separations of 6 - 435 AU. 



3.7. Praesepe 

The Praesepe open clu ster is located »T70 pc from the 
Sun dHamblvetalJll995h and has an age of a bout 600 Myr 
dBoudreault et all 1201 oT IKraus & Hillenbrand! d2007l) iden- 
tify 1050 ± 30 Praesepe members earlier than spectral type 
M5, corresponding to 550 + 4OM . The mass-function in 
Praesepe agrees well with the present-day MF for field stars 
(IKraus & Hillenbra ndll2007l: iBoudreault etaD2010t iBaker et al.1 
l2010h . lAdams et al.ld2002l) state the half-mass radius to be 1 .°25, 
corresponding to a; 4 pc at the distance of Praesepe, and a total 
ma ss of » 6QQM within 4° = 11.8 pc. 

IPatience et al. (2002) find 12 binaries among 100 Praesepe 
members in their sample with separations of » 1 8 - 390 AU. The 
binar y fraction is thus 12 percent, about half as large as in the 
field. lBouvier et aT (2001), in contrast, derive a corrected binary 
frequency of 25.3 ±5.4 percent, which is in quite good agreement 
with that measured in the Galactic field. Their work, containing 
26 binaries between 15 and 600 AU among their larger sample 
of 149 cluster members, is used here. 



3.8. Pleiades 

The Pleiades open cluster is located at a distance of about 

1 L 1 I _ 1 

125 pc (Robichon et al. 1999) and ha s an age of about 70 - 
120 Myr dStaufferet al.l ll998; M artin et al.lll998l) . Its approxi- 
mately 1245 member stars have a total mass of about 74OM 
dConverse & Stahlerl2008"HPinfield et all 19981) and its half-mass 
radius is 2 pc dRaboud & Mermilliod 1998). Dynamical mod- 
els of the Pleiades show that the IMF of the Pleiades might 
have been similar to the Galactic field-star (i.e. canon ical) IMF 
for stellar masses below 2M dMoraux et al.1 120041) and that 
it may have a precursor similar to the ONC ( Kroupa et al.l 
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1200 ll). Its initial binary fr action might have been close to unity 
(IConverse & Stahlerll20ip|) and its initial dis tribution of periods 
si milar to that in T a urus dKroup a et al. 2 001b . 

iBouvier et al.1 (1 19971) studied the binary properties of 
144 systems in the Pleiades finding 22 binaries and 3 triples 
in the orbital separation range 10 - 900 AU. The corrected 
binary-fraction of 28 + 4 per cent agrees well with the Galactic 
field value of 27 per cent over the same range. 



4. Results 

4. 1 . The young cluster sample 

Figs. [2] E] and [4] show the results of our analysis. Each best-fit 
theoretical distribution was evolved for the time indicated in the 
upper right comer of each panel. This time roughly corresponds 
to the age of the region (Sect. [3]l or if the age is uncertain or 
shows a spread it is an estimate of the upper age limit. 

Although the seven regions exhibit present-day stellar den- 
sities ranging from 1-10 stars pc~ 3 (Taurus) to 2 x 10 4 stars 
pc~ 3 (ONC), the model is able to reproduce the observational 
data for all of them (though only poorly for USco-A, although 
see below), requiring solely different initial densities. All dynam- 
ically evolved initial binary populations are likely parent func- 
tions of the respective observational data according to the % 2 - 
test. 

It is found that the more depleted the binary population with 
respect to the initial distribution, the higher the required initial 
stellar density, p ec \, to reproduce the observational data (for the 
same age of the population). Taurus (Fig. [2] upper panels) is 
very sparsely populated and still contains a large binary popu- 
lation. Therefore, the initial stellar density was also low, to en- 
sure that the binary disruption was inefficient. According to the 
model, Taurus sub-clusters formed with p e ci = 35OM pc~ 3 , but 
is consistent with a semi-major axis distribution close to the ini- 
tial distribution within the confidence intervals (Fig. [2] upper left 
panel). The current masses and sizes of the individual clusters in 
Taurus are indeed comparable to the inferred value of p ec i (Fig. [2] 
upper right panel). 

In contrast, the similarly aged p Ophiuchus (Fig. [2] lower 
panels) and ONC (Fig. [3] middle panels) show a much stronger 
depletion of binaries and therefore require higher birth densi- 
ties for more efficient stimulated evolution. The present-day val- 
ues are thus farther away from the initial conditions constrained 
here. This is readily understood as a result of the stronger expan- 
sion within the same time of initially denser regions owing to the 
more effective two -body relaxation and bin ary burning (fig. 4 i n 
Marks et al. 201 1) and/or expulsion of gas ( Kroupa et al.ll2Q0lT) . 

According to our model, the binary fraction above 1000 AU 
is two percent for the ONC. For a population of cur- 
rently approximately 3500 systems, this would correspond to 
many more binaries in this semi-major range than t he up- 
per l imit of three binarie s estim ated by Scallv etail |l999) 
for a lDuquennov & Mavorl dl99ll) -tvpe distribution. The model 
indeed suggests a semi-major axis distribution that is more 
strongly depleted in binaries (Fig. [3] middle panels). A detailed 
study including instrumental limitations an d projection effects 
would need to be performed to address the IScallv et"aT1 (1 19991) 
analysis. This is however beyond the scope of the present study. 

Our model provides only a poor fit to the data for USco A. 
Although a fit formally exists, the large confidence intervals 
would allow a large range of values for the initial density of 
USco A. It might thus be difficult to explain this data in terms 
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Fig. 5. Comparison of present-day and inferred initial stellar 
densities of typical regions. All seven regions evolve to lower 
densities away from their initial values. Ranges for the present- 
day densities are calculated from the total stellar masses and 
half-mass radii reported in Sect. [3] USco-A is not plotted here as 
no present-day radius with a given enclosed mass is reported in 
the literature that could reasonably be compared to the half -mass 
densities of its initial configuration constrained here. Note that 
some mass estimates only cover a limited spectral-type range. 
The filled circl es indicate the deduced initial conditions for some 
Galactic GCs dMarks & Kroupal2010l) . 



if dynamical modification of the Kroupa ( 1995b) initial distribu- 
tion. USco A is a complicated region, part of an OB association, 
and already far more extended than a typical open cluster. It is 
likely that USco A may contain projected members of the asso- 
ciation. The USco-B subgroup in the Scorpius-Centaurus asso- 
ciation is not analys ed here, as the binary fra ction increases with 
increasing period dBrandner & Koehlerll998l) . This is incompat- 
ible with the dynamical evolutio n of its binary population an d 
may be a contamination effect dKroupa & PeuMj otzens 201 lb . 
We note, however, that iParker & Goodwinl d2012l) caution that 
binary disruption is a highly stochastic process that influences 
the appearance of observed binary distributions. 

The data for the open clusters Praesepe and the Pleiades 
(Fig-SJ 0.6 and 0.125 Gyr respectively) are also closely repro- 
duced by the model0 The initial distribution for the best-fit p ec \ 
was evolved for only 5 Myr, but by then the evolutio n of the bi- 
nary distributions had effectively come to an end (Marks et al. 
1201 ll) so that a comparison with the much older regions should 
be possible. That the binary population does not evolve signif- 
icantly after 3-5 M yr is clearly evident from a comparison of 
iMarks et all d201 lb with lKroupa et alj d2001l) . 

According to the model, for IC 348 and the Pleiades the best- 
fit initial densities are very similar, as are those for Praesepe, 
p Ophiuchus, and Chamaeleon (Fig. [5] not considering the rather 
large confidence intervals). Taurus and the ONC occupy the low- 
and high-density end, respectively. The large confidence inter- 



2 Note that larger values of P{> xt&o) ^° not signify a statistically 
tighter fit than for the younger regions, but that one would be less con- 
fident that the theoretical distribution is not a parent function were the 
null-hypothesis to be rejected. 
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vals even allow Praesepe and the Pleiades to have formed with 
the ONC density, as discovered before by iKroupa et al.l (1200 ll) 
using independent arguments. 

4.2. Cluster evolution depending on density and comparison 
with the initial conditions of globular clusters 

Mark s & Kroupal {2010) provide independently obtained con- 
straints of the initial conditions for some Galactic globular clus- 
ters (GCs). They use the imprint of residual-gas expulsion still 
visible in the present-day stellar mass function of GCs to back- 
trace the present cluster properties to their initial values. Their 
results neatly join the densest young clusters investigated in 
Fig. [5] suggesting a connection between the initial and present- 
day densities of clusters, althoug h this is difficult t o infer from 
the young cluster sample alone. Mark s et al.l (1201 ll) indeed find 
that a dense initial configuration can be less dense today than 
another initially less dense cluster. 

Fig. [6] shows that the constraints compare well with masses 
and d ensities of mo l ecular clumps that have just begun to form 
stars iMueller et al.l (|2002|) : IShirlev et ail (120031) : iFontani et all 
J2005ir The masses of the clumps from these sources have been 
multiplied with a star-formation efficiency of one-third for com- 
parison with the stellar masses and densities calculated here. 

Except for Taurus, all of the investigated regions lie close 
to or above the threshold of the sur face density of massive-star 
formation of 1 g cirT 2 sugg ested by Krum holz & McKed ( I2008L 
which is the same for all cluster masses). This value corresponds 
to 4788M pc~ 2 and, assuming an average pre-cluster molecular 
cloud diameter of 0.76 pc, the average clump size in the three 
aforementioned observational data sets, corresponds to p ec i = 
4788/0.76 = 6202M G pc~ 3 (thin solid line, for a homo geneous 
mass distribution within the clumps). 

According to this interpretation, far more-massive and 
denser clumps were needed to form proto-GCs. Considering our 
derived constraints (solid squares in Fig. [6} and the data sets 
for GCs (solid circles) we find a relation between stellar mass- 
density and stellar mass of the form 

logioPeci = a x log 10 M ed + b , (6) 

where a = 0.61 ± 0.13 and b = 2.08 ± 0.69 (thick solid line). 
Fitting only to the young cluster data a = 0.59 ± 0.22 and 
b = 2.26 ± 0.57 (thick dashed line), i.e. both fits are indis- 
tinguishable. We note however that the trend without the GCs 
is mainly driven by Taurus and the ONC. Excluding these two 
only a rather fiat initial stellar mass-density-stellar mass relation 
emerges. Equation [6] implies that there is a weak stellar mass- 
half-mass radius dependence of star clusters at birth of the form 



£h 
pc 



3(M ecl /M B ) 1 
8tt ■ 10* 



1 +0 07 x 

U - iU -0.04 x 



M ec i 



U.13±0.04 



(7) 



which is consistent with a relation between the effective ra- 
dius and star cluster mass for a sample of clusters i n non- 
interact ing spiral g alaxies, R e ff oc mQ- 1q± Q3^ found bv l Larsenl 
(120041) . IZepf et alJ d 19991) and IScheepmaker et al.l (|2007|) simi- 
larly demonstrated that there is a mass - luminosity relation for 
clusters that is shallow for both the galaxy merger NGC 3256 
and young clusters in M5 1 . 

In Fig. [71 the clusters that we have studied appear to follow 
an evolutionary sequence that depends on the number of com- 
pleted initial crossing-times, 



'crjnitial 



Vg 



eel n 



(8) 



The Galactic GCs dMarks & Kroupal feOlO) are older and, thus, 
dynamically more evolved but show similar initial-to-present- 
day density ratios as Chamaeleon, Praesepe, and the Pleiades. 
The molecular clumps presented in Fig. [6] are expected to lie 
at the bottom-left corner of Fig. [7] as these continue to collapse 
before star formation sets in. 

The results presented here suggest that bound star clusters 
reduce their birth density by a factor of a few hundred within 
10 - 100 initial crossing times and then evolve at about constant 
density. The standard secular evolution trac ks (without gas and 
binaries) as predicted by A^-body models ( Kupper et al. 2010h 
with initial conditions as stated in the caption of Fig. Q are un- 
able to explain the constrained location for the young clusters 
and GCs. The energy that initially expands clusters consider- 
ably instead needs to be produced to achieve consistency with 
the young cluster and GC data presented here. This production 
may be possible by residual-gas expulsion and/or the burning 
of binaries, which enables the freshly-hatched clusters to as- 
cend the track described by the real clusters in this work. The 
most massive clusters may evolve into old GCs, but, after their 
gas expulsion- and binary -burning-driven expansion, at a slower 
pace in terms of the evolution of their density, which is probably 
similar to the quiescent evolution of the presented A^-body mod- 
els. In this sense, GCs and open clusters (Pleiades, Praesepe) 
may have formed v ia a single clu ster formation mode, as previ- 
ously suggested bv lLarsenl d2002l) . The intrinsic slope evident for 
the GCs in Fig.|7]may however suggest that the picture is more 
complicated. This needs to be tested by self-consistent A^-body 
models of star clusters including gas expulsion and binaries as 
cluster heating sources. 

5. Discussion 

We have shown and discussed that evolving an assumed univer- 
sal initial binary population in star clusters for the age of the con- 
sidered object leads to theoretical semi-major axis distributions 
that are consistent with the observed ones in six young clusters 
and star-formation regions, as well as in two older open clusters. 
Only the stellar densities at birth we needed to carefully select 
to adjust the strength of stimulated evolution. For all investigated 
objects, the theoretical semi-major axis distributions turn out to 
be p arent functions of observational data. 

Mark s et aD (1201 ll) demonstrated that there is an age-density 
degeneracy, i.e. regions of different age can have the same bi- 
nary population if their birth densities are properly adjusted. 
This degeneracy is lifted here because the investigated popu- 
lations have a reasonably well-known age allowing IDPS to to 
be applied to each cluster so that the initial density of a typical 
regions' progenitor (i.e. a » 0.5 pc sub-clump or individual clus- 
ter in the case of Taurus but the whole ONC) can be inferred. 
Taurus appears to host a largely unevolved bi nary population, in 
agree ment with explicit A^-body modelling dKroupa & Bouvieil 
2003). The present-day observed masses and sizes of individual 
~ 0.5 pc clumps in Taurus are therefore close to what is expected 
based on IDPS of the regions' initial conditions. In the case of 
the ONC we need to assume the highest density at birth to turn 
the initial separation distribution into the observed one and its 
present-day mass and radius are rather far away from the ones 
at birth. In-between these two extremes, IC 348 and the Pleiades 
on the one hand and Praesepe, p Ophiuchus, and Chamaeleon 
on the other share very similar initial densities according to this 
analysis. For each cluster, the calculated initial density agrees 
with the pre-stellar cloud clump densities observed in molecu- 
lar clouds (Fig. |6j. This therefore suggests that the expansion 
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caused both by residual-gas expulsion and efficient energy gen- 
eration owing to binary-burning increases with the ratio of age 
to t a 

.initial 

(Fig.0. 

Furthermore, by adding th e constraints on the birth densities 
for Galactic GCs obtained by iMarks & Kroupal d2010l) . we un- 
cover a stellar mass - half-mass radius relation of star-forming 
cloud cores (eq. |7). This relation implies that the radii of birth 
clusters have a weak stellar mass dependence, r^ oc M° - 3±0 04 , 
that is consistent with other observation al data (IZepf et al. 1999; 
Larsen 2004; Scheepmaker et al . 20071). 

iKing et alJ d2012h compare the binary properties of five 
young clusters analysed here (except USco-A). They find no cor- 
relation between the present-day binary fraction and the present- 
day (instead of initial) number density, at least if Taurus is not 
considered. The present-day density might however not neces- 
sarily reflect the initial conditions as the large density of an ini- 
tial configuration (with strong binary burning) can drop below 
the density of an init ially less dense co nfiguration within only a 
few Myr (cf. fig. 4 in lMarks et al.ll20 lib . Residual-gas expulsion 
may enhance this effect. 

Using A^-body computations, IKing et all d2012h tentatively 
argue for fractal initial conditions in these regions, even though 
this is not the best A^-body model for all their studied regions. 
They invoke this substructure on different scales (< 1 pc for 
ONC- and IC348-type clusters, but < 10 pc for Taurus-type con- 
figurations). The sizes of the clumps in Taurus are however com- 
parable to the extent of the ONC and IC348, as evident in their 
figs. 7-11. Thus, there appears to be one dynamically relevant 
scale, i.e. the scale on which dynamical (binary) evolution is im- 
portant, of about w 0.5 - 1 pc. This size corresponds well t o the 
sizes in the A^-body models on which the Marks et al. (201 1) an- 
alytical model is based and that are used here to derive the initial 
densities. 

The analysis here finally suggests that stimulated evolution 
plays an important role in shaping the present-day properties of 
binary populations in star clusters. There exists of a formal initial 
binary distribution function that is consistent with the available 
multiplicity data of pre-main sequence and Class I protostellar 
populations ( Krou pa & Petr-G otzens 2011) and that helps us to 
understand the binary distributions in the eight objects investi- 
gated here simultaneously, assuming dynamics only. This sug- 
gests that there are no significant variations in binary-star for- 
mation depending on the ambient condition required, but that 
variations in the binary populations are mainly caused by stellar 
and dynamical evolution. This agrees with previous work and 
such an assumption wou ld be similar to the universality hypoth- 
esis o f the stellar IMF (iKroupal I20HI iKroupa & Petr-Gotzensl 
l201lh . 

The alternative is that the birth binary populations depends 
on the cloud density such that a larger fraction of binaries form 
in low-density environments. This is a trivial solution to the ob- 
served trend. However, the IMF is k nown to be invariant for the 
star-forming conditions probed here (Kroupa et al. 201 1), and it 
is unclear how the IMF of systems could be separated from the 
IMF of all individual stars if there were a variation in the prop- 
erties of the initial binary population. 
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Fig. 2. Left column: Semi-major axis distribution functions. Both panels show the initial dKroupal 199 5b) semi- major axis distribu- 
tion i n star clusters (dashed curve), as well as the canonical solar-type Galactic field distribution (solid curve, Duquennov & Mavor 
11991b . The black squares are observational data from the references given in the legend to which the model is fitted (in Figs. [2]|4] 
corresponding to the most recent data or the largest sample). Vertical errorbars indicate y/N errors, the horizontal ones indicate 
the width of the bin. The initial semi-major axis distribution evolves into the solid histogram for an initial stellar density within the 
half-mass radius, p e cb as indicated in the legends. The respective ^ 2 -tests yield the indicated X^m' anc ' P(- ^min-*" va ^ ues ' ^ tne initial 
distribution is evolved for the time t in a cluster of density p ec i (note that t is mostly an upper limit to the age of the population). 
The resulting distributions are likely parent functions for the observational data. The dashed histograms are 10 percent confidence 
limits, outside which a distribution would be rejected as a probable parent function with > 90 percent confidence. The overall model 
binary-fraction, /bin (integrated over the whole semi-major axis range), is indicated in the legend. Right column: The thick solid 
and dashed lines are lines of constant volume density within the half-mass radius according to the best-fit p ec \ and the confidence 
intervals, respectively. Thus, the model suggests that the total stellar mass and half-mass radius in a typical, dynamically active 
region (e.g. a sub-clump) of the star-forming environment's progenitor lies somewhere in the grey shaded region. For comparison, 
observed present-day masses and spatial extensions are shown (Sect. [3]). Note that the masses and sizes shown denote different 
spatial scales (e.g. the constant density lines depict the typical total stellar masses, M ec i, and the half-mass radius, r n , required to 
produce the present binary population, while the symbols might refer to different quantities as indicated in the legend - i.e. the total 
stellar+gas mass within the whole region, not solely individual clumps). 
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Fig. 4. As Figs. [2] and [3] but for USco-A, Praesepe, and the Pleiades. Although Praesepe and the Pleiades are much older (0.6 and 
0.125 Gyr, respectively), the initial distribution is evolved for t — 5 Myr only when binary-burning has largely ended (see the text). 
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Fig. 6. Constraints on the initial volume-densities within the half-mass radius and masses derived in this work for the seven clusters 
(filled squares) versus the initial stellar mass. The indicated errors in mass correspond to the observationally inferred present-day 
mass on the left end of a bar and two times the present-day mass on its right end, to be understood as an estimator of the possible 
initial-mass range. F illed circl es are Galact i c GCs as in Fig. [5] Underlaid as grey sy mbols are data of mole cular cloud clumps of 
iMueller et al. (2002, crosses). IShirlev et al.l (l2003[ squares), and lFontani et alJ (I20051 circles) as collated by Parme ntier & Kr oupa 
(EoT ]]). These are known to have already begun forming stars. The clump masses have been multiplied with a star-formation 
efficiency of one-third to compare to the stellar mas ses and densities in f erred in the present work. The thin solid black-line is the 
threshold for massive-star formation evaluated bv iKrumholz & Mc Kee (2008, 1 g cm 2 , see the text). The thick solid black-line 
is a least squares fit to both the young cluster and GC data (eq. implying that there is a mass-radius relation (eq. [7} for star 
cluster-forming cloud clumps. The dashed line shows the result when the GCs are excluded from the fit. 
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Fig. 7. The ratio of the constrained initial to the present-day density as a function of the number of completed initial crossing-times, 
fcr,initiab for seven clusters (filled squares, where USco-A is not shown for the same reason as it is not depicted in Fig. |5). Errors in 
the x- and y-directions follow from the error in p ec i and error propagation. These clusters may originate from the bottom-left corner 
of the diagram, the region supposedly occupied by the molecular clumps shown in Fig. [6] The clumps are probably still collapsing 
and have not yet reached their peak d ensities, hence have a dens ity ratio < 1 . Filled circles indicate the constraints for Galactic GCs 
as in Figs. [5] and [6] (ages taken from iMarin-Franch et al . 2009). The solid arrows indicate a possible evolutionary sequence. This 
evolution cannot be understood in terms of single-star A^-body star cluster models without gas exp ulsion. These models are shown 
by the evolution of the M ec \ =2x1 4 M clusters with the indicated initial-^ and a | Kroupal d200ll) IMF bet ween 0. 1 M and 1 .2M 
on circular orbits around a realistic Galactic potential over 4 Gyr (Kiipp er et al. i i2oior solid black curves). Residual-gas expulsion 
and/or energy generation through binary -burning probably needs to be invoked to explain the expansion of the real clusters (see the 
text). 
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